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Cosmic Acceleration from Elementary Interactions 
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It is possible to generate an accelerated period of expansion from reasonable potentials acting 
between the universe particle constituents. The pressure of primordial nucleons interacting via a 
simple nuclear potential is obtained via Mayer’s cluster expansion technique. The attractive part 
of the potential engenders a negative pressure and may therefore be responsible for the cosmic 
acceleration. 
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I. INTRODUCTION 

Cosmic primordial acceleration is usually derived from 
a scalar field. The idea of using a scalar field arose in the 
context of GUT’s |jj, at first with the Higgs boson field 
in mind. This attribution has been dropped in the long 
run, as the Higgs field did not reproduce all necessary 
features required by the inflationary period fgj. As a 
consequence the inflaton - an ad hoc scalar field leading 
to the desired properties - was introduced. This arbitrary 
field, however, has no simple interpretation in terms of 
fundamental physical phenomena — which justifies the 
search for alternative explanations. 

We intend here to show that primeval acceleration may 
come from the strong interactions in the constituents’ 
equation of state. Such short range interactions are taken 
into account via the well-known and systematic method 
of cluster expansions j 3 |. 

An application to cosmology of a toy-model equation 
of state (EOS) for systems with interaction has been 
made in jj] , where the elementary constituents of the uni¬ 
verse were taken to be hard-spheres. A recent paper |j 
has proposed an explanation for the inflationary period 
through a short range interaction, just what we shall do. 
Nevertheless, the modification on the equation for the 
pressure was not obtained through the Mayer method 
here employed, which brings to light a time scale deter¬ 
mined by the temperature and leads to a natural passage 
from the accelerated period to the decelerated one. 

In section HD we synthesize the main characteristics of 
our model, specifying the relevant particles of the matter 
content. Section mu deals with the construction of the 
equation of state for the interacting media using Mayer’s 
approach of virial cluster expansion. The nuclear inter¬ 
action is modeled by the effective potentials presented in 
section UYI The results for our simplistic model point to 
a period of accelerated expansion and are discussed in 
section El 
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II. ON THE PRE-NUCLEOSYNTHESIS PERIOD 

The pre-nucleosynthesis period is the stage immedi¬ 
ately preceding the cosmological formation of the light 
elements (deuteron, He +2 , etc). This formation only 
started when the universe mean energy attained values 
around the deuteron binding energy ( Eb — 2.23 MeV). 
So, we roughly characterize this period by the red-shift 
z>2x 10 12 or, in terms of energy, kT > 4 MeV. 

Radiation dominates the universe content in this 
phase. Photons generate a panoply of particles via pair 
production, with populations dependent on the energy 
of the thermalized system HQ In particular, at hun¬ 
dreds of MeV nucleons produced by this process are in 
chemical equilibrium Q. At these energies the numeri¬ 
cal density tin of nucleons deriving from the radiation is 
much superior to that (rib) of the present-day nucleons, 
which can be consequently neglected. 

A realistic scenario takes into account, besides the nu¬ 
cleons, light particles appearing even before the nucle¬ 
ons in the pair production process. Examples are pions, 
eletron-positron pairs, muons. For simplicity, we will 
adopt as relevant content photons, nucleons and anti¬ 
nucleons (protons and anti-protons, neutrons and anti¬ 
neutrons). We will neglect eletrons and positrons, pions, 
neutrinos, muons and anti-muons. 

The presence of charged particles would require the 
consideration of the electromagnetic interaction among 
the constituents as a source for the gravitational field. 
However, tin » rib implies symmetry between matter 
and anti-matter. We shall suppose overall Debye screen¬ 
ing: in large enough scales, the system will be electri¬ 
cally neutral. Once the electromagnetic interaction is 
neglected, we are left with the strong and weak interac¬ 
tions. Obviously, the weak nuclear interaction is much 
less intense compared with the strong one: the strength 
of the former is 10" 13 orders of magnitude of the last 
[H. This justifies the approximation we will adopt: our 
equation of state for the pre-nucleosynthesis period con¬ 
siders only the strong nuclear interaction. In short, our 
model (7 + NN) is constituted by photons plus interact¬ 
ing nucleons treated as classical non-relativistic particles 
obeying Boltzmann statistics. 
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III. EOS OF A SYSTEM UNDER 
INTERACTION 

With the hypothesis introduced above, the equation of 
state of the system composed by photons ( 7 ) and nucle- 
ons/antinucleons (N/N) will have the form: 

p(kT, n N ) = p^(kT) + p N (kT,n N ), (1) 

where kT is the energy and n n is the 
nucleon+antinucleon numerical density. 

We shall approach the calculation of the strong inter¬ 
action between nucleons by Mayer’s method Q. In this 
formalism, the pressure of a classic gas of interacting par¬ 
ticles is calculated through the virial expansion: 

p N (kT,n N ) = n N kT y ^ai(T) ( nN ^ ^ , (2) 

1=1 \ 9 J 

cii(T) being the virial coefficients, viz. 


terms of this parameter. This requires the determina¬ 
tion of tin = niy(kT). Supposing equality between the 
nucleon and antinucleon concentrations and masses, the 
numerical density will be 
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where pn is the nucleon chemical potential. 

The system in equilibrium must assure the formation 
of the nucleons, i.e., p^ = —m^rc 2 . 


IV. GRAPHICS OF p(kT) FOR 
PHENOMENOLOGICAL POTENTIALS 


In a simplistic manner, we describe the strong interac¬ 
tion by a nuclear potential of the square-well+hard-core 
type shown in Figure ^ The choice of this potential is 


ar = 1; a 2 = - J J (e~ u ^ kT - 1 )r 2 dr; (3a) 

0 

00 

03 = 3A ^ JJ ^ 12 ^ 13 ^ 23d3ri2rf3ri3; ( 3b ) 
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and so on, given in terms of the mean thermal wave¬ 
length A = h/yJ^-nm^kT of the nucleons, and of the 
two-particle Mayer function fy = e _/3,ii5 ' — 1 which de¬ 
pends on f3 = 1 /kT and on the inter-particle potential 
Uij = u(\fi — fj |). The second virial coefficient a 2 in <!3bl) 
reflects the interaction of the particles two by two (the 2 - 
cluster); and 03 accounts for the interaction among three 
particles, g is the number of possible values for the in¬ 
ternal degrees of freedom. Four species are to be taken 
into account (protons, anti-proton, neutron and anti¬ 
nucleon). Besides, our species have spin A: g = 4x2 = 8 . 

Under our assumptions, only short-range nuclear forces 
will be at work. It is consequently reasonable to admit 
that interaction occurs only between few particles, and 
that pw as given by © can be conveniently approximated 
by its first three terms. This hypothesis is reasonable as 
long as the density is not too high. Substituting m in 
the equality © leads to the EOS 


p(n N , kT) = 


45 h 3 c 3 


(kT) 4 + n N kT+ 


+riNkT 


a2{kT) (^] +a3{kT) (^\ 


V 9 


9 J 


( 4 ) 


The first term, p 1 (kT ), comes from the black-body ther¬ 
modynamics. 

The natural time parameter along the thermalized pe¬ 
riod predating recombination is just kT. To examine the 
history of Eq. 0 , it would be convenient to express it in 




FIG. 1: Scheme of the supposed strong interaction potential. 


justified by two facts: (i) it shows a behavior qualitatively 
similar to some successful phenomenological nuclear po¬ 
tentials, e.g., those discussed in EULGl ; (ii) it makes it 
possible to calculate analytically the corrections to the 
ideal case. 

The potential parameters are determined from the 
deuteron binding energy, from its mean square radius — 
which constrain the well width (b — c) = 1.3 fm and its 
depth Vo = 75.6 MeV — and from nucleons scattering 
data at high energies — which determinate the hard core 
extension c = 0.4 fm (cf. reference @). 

Function p(kT) is obtained by substituting © in © 
and using the a 2 (fcT) and a^(kT) explicit expressions. 
Its aspect is shown in Figure |5] It exhibits an increas¬ 
ing behavior at low energies (kT < 330 MeV). This is 
only natural, since at these energies the nucleons density 
is still too small to cause a relevant interaction effect. A 
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FIG. 2: Pressure of the interacting nucleons plus photons as 
a function of the energy: given in units of MeV/ A 3 , where 
Ac = ftc/mjvc 2 , and the energy is given in MeV. 


FIG. 3: Plot of p + 3p as a function of kT in MeV. In our 
units, AnG/3c 2 = 1. Recall that the cosmic time and the scale 
factor decrease as the temperature increases. 


sudden change occurs, however, at around 375 MeV. The 
pressure starts decreasing as the energy rises, becoming 
negative at around 420 MeV. This surprising behavior is 
due to the action of the interaction term d 2 (kT) which 
dominates in ©• To be more specific: it is the attrac¬ 
tive part (negative sector) of the potential in a 2 (kT) that 
surpasses all the other terms. 

It is worth noticing that for energies until 490 MeV 
the two by two interaction term (a 2 (fcT)nivA 3 ) between 
nucleons is at least five times larger than the three by 
three interaction term (ci^(kT)n 2 N \ Q ). This is a point in 
favor of truncating the series 0 in the first terms. 


Energy (MeV) 

Thermal history phase 

0 < kT < 100 

decelerated expansion (p = f) 

100 < kT < 450 

transition phase (— % < p < 

450 < kT < 470 

accelerated (—p<p< — f) 

kT > 470 

ghostly accelerated (p < —p) 


TABLE: Thermal history phases of the model with inter¬ 
action. Notice how the model links naturally the accel¬ 
erated expansion period with a radiation type decelerated 
era. 


V. FINAL REMARKS 

As the pressure, the energy density pw(kT,nN ) may 
be expressed in terms of a series of type 0 . The method 
is the same presented in section m and the result is: 

p(kT,n N ) = 15 ^ 3 e 3 (kT) 4 + n N ^m N c 2 + —^j + 

T <9a,3 \ njvA 3 
a3 ~ HItt) g \ ’ 
( 6 ) 
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ton being, we recall, the nucleon mass. 

With the pressure equation 0 and the equality 0 , 
we can rewrite the Friedmann equation 
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in terms of kT. The behavior of p + 3p is given in Figure 

01 

Combined with Q, it shows that positive acceleration 
occurs from 450 MeV on. A detailed analysis of the 
Friedmann equations brings forth a division in the ther¬ 
mal history of the model with interactions in four distinct 
periods: 


Further improvement would come from taking into ac¬ 
count the neglected particles. The nuclear interaction 
could also be approached in a more sophisticated way, 
e.g., by using the phenomenological potentials given in 
refs. El introducing quantum corrections a la 

Bethe-Uhlenbeck ||, adding relativistic corrections, etc. 

The idea that an attractive potential between the con¬ 
stituents can led to accelerated expansion can be tenta¬ 
tively applied to the present day observed acceleration. 
The scenario is then that of matter interacting through 
the gravitational potential. The latter being attractive, 
the expected global result would be a decreasing in the 
value of the pressure towards negative values, eventually 
causing the acceleration. Around the same period, but 
locally, the gravitational potential would became suffi¬ 
ciently effective to engender the large structures in the 
universe. This is consistent with the observational data, 
which indicate that both the accelerated expansion E3 
and the large scale structure formation 0 ] have begun 
recently. 
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